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We present an overview of the measured transport properties of the two dimensional electron 
fluids in high mobility semiconductor devices with low electron densities, and of some of the 
theories that have been proposed to account for them. Many features of the observations are not 
easily reconciled with a description based on the well understood physics of weakly interacting 
quasiparticlcs in a disordered medium. Rather, they reflect new physics associated with strong 
correlation effects, which warrant further study. 
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As low density two dimensional (2D) electronic sys- 
tems with increasingly high mobility have become avail- 
able, there has accumulated experimental evidence of a 
set of low temperature phenomena which cannot be eas- 
ily understood on the basis of traditional (Fermi liquid 
based) metal physics. 1 More precisely, these are systems 
with a large ratio between the typical potential and the 
kinetic energies, r s = 1 / \j 'im{a-Q) 2 , where n is the areal 
density of electrons and = fi 2 e/m*e 2 is the effective 
Bohr radius. A generic feature of electronic systems is 
a low temperature evolution as a function of n from a 
conducting "Fermi liquid" (FL) state at large n to an 
insulating "Wigner crystalline" (WC) state at low. To 
obtain as clear a perspective as possible on this physics, 
we focus on 2D electron (and hole) systems in which r s is 
large, so interaction effects manifestly cannot be treated 
perturbatively, but in which the system remains a fluid, 
hence thermodynamically distinct from a WC. 

It should be stressed that while the experimental data 
exhibits striking features which have been confirmed by 
multiple groups on various materials, there is disagree- 
ment concerning the correct theoretical perspective from 
which to view these experiments. Of the unsettled is- 
sues, the most vexed is the proper theoretical interpreta- 
tion of the experimentally observed metal-insulator tran- 
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sition. There is also much that remains to be understood 
concerning the observed transport anomalies in samples 
which arc far from the critical region. Here we summa- 
rize the salient experimental facts, stressing the anoma- 
lous character of the results obtained in conditions both 
close to and far from the metal-insulator transition. We 
also summarize some of the major theoretical approaches, 
including a discussion of both their successes and short- 
comings. The primary purpose of this Colloquium is to 
stimulate further experimental work on this subject, par- 
ticularly work aimed at unraveling the non-critical be- 
havior of the metallic and insulating states that occur 
in strongly correlated two dimensional electron fluids far 
from the metal-insulator transition. 



I. EXPERIMENTAL SIGNATURES OF NON-FERMI 
LIQUID BEHAVIOR 



In this section, we summarize some of the experi- 
mental results on large r s 2D electron and hole gases 
(2DEG's and 2DHG's) in highly conducting semicon- 
ductor heterostructures that we contend are incom- 
patible with a Fermi liquid based theory and with 
the single particle theory of weak localization. We 
will particularly focus on experiments on the follow- 
ing systems: n-Si MOSFETs, 2 p-GaAs hetero-junctions 
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of the critical value of r s at which the energies of a uni- 
form FL and WC are equal is r* ss 38, while for the 
devices we have in mind, r s ~ 5 — 20 in the case of the 
Si MOSFETs and ~ 10 - 40 for the p-GaAs devices. 



One of the key points to notice as we discuss the data 
below is the great similarities in the structure of the data 
from the different devices. This implies that the observed 
anomalies represent robust "universal" behaviors of the 
2DEG, which are largely independent of details. This is 
correct in spite of the fact that there are significant dif- 
ferences between the electronic structures of the various 
devices. In p-GaAs and n-GaAs heterostructures, the 
electrons or holes occupy a single band with an isotropic 
effective mass, while in n-Si MOSFETs, there are two de- 
generate valleys and correspondingly a non-trivial struc- 
ture to the effective mass tensor. In the heterostruc- 
tures, the interactions between electrons at large separa- 
tion are Coulombic, while in MOSFETs, the interaction 
between electrons is dipolar at distances large compared 
to the distance to the metal gate. In p-GaAs and p- 
SiGe, the spin-orbit coupling may be significant, while in 
n-GaAs and Si MOSFETs it is clearly insignificant. In 
Si MOSFETs the disorder potential is believed to have 
short-ranged correlations, while in most of the other 2D 
systems considered here, it is believed to be long-range 
correlated (i.e. due, primarily, to distant charged impu- 
rities). 
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A. Near critical samples with p ~ h/e 2 

Both the 2DEG and the 2DHG exhibit a zero- 
temperature metal-insulator transition as a function of 
the electron density n. This effect has been observed 
in all materials in which large r s samples can be made 
with sufficiently high mobilities: p- and n-Si MOSFETs, 
p- and n-GaAs quantum wells, p- and n- SiGe quantum 
wells, and n-AlAs quantum wells. Examples of the ex- 
perimental data showing the resistivity p(T) for different 
electron concentrations in various systems are presented 
in FigED 

In all cases, at low temperatures T, the resistiv- 
ity p(T) exhibits a "metallic" temperature dependence 
(dp(T)/dT > 0) for electron concentrations n in excess 
of a well defined critical value, n c , and dielectric behav- 
ior (dp(T)/dT < 0) for n < n c . Moreover, typically 
at the lowest temperatures p <C h/e 2 for n > n c , while 
p ^> h/e 2 for n < n c . (Note that, were the Drude formula 
valid here, p = h/e 2 would correspond to a mean-free- 
path, £, equal to the Fermi wave-length: kp£ = 2tt. The 
quantum of resistance, h/e 2 , is thus the upper limit for 
the possible regime of applicability of Boltzmann trans- 
port theory) 

An issue that has been debated is whether there is an 
actual transition or just a rapid crossover. This issue 
is difficult to resolve unambiguously since a continuous 
metal-insulator transition is sharply defined only at zero 
temperature as the point at which the resistance changes 
non-analytically from being finite (metallic) to infinite 
(insulating). (Note, under the heading "metallic" we 
include the possibility of a "perfectly metallic" phase 
in which p is non-zero at any non-zero T, but p — > 
as T — > 0.) However, at non-zero temperature, even 
where there is a zero temperature transition, at any 
T > 0, such a transition would manifest itself as a rapid, 
but analytic crossover. It is clear from FigQ] that the 
resistivity changes by many orders of magnitude as n 
varies over a modest range near n c . The lower the 
temperature is, the more violent the resistivity changes. 
For comparison, in FigfT] (f ) we present data for the 3D 
metal-insulator transition in Si:P. No one doubts that 
this data reflects a zero temperature metal insulator 
transition, although from a purely empirical viewpoint, 
the evidence is no better (and no worse) than in 2D. 
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FIG. 1 Critical behavior of the resistivity near the 2D metal- insulator transition in a Si MOSFET (a), p-GaAs/AlGaAs 
heterostructure (b), n-SiGe heterostructure (c), n-GaAs/AlGaAs heterostructure (d), and p-SiGe quantum well (e). The 
corresponding ranges of r s are 14 — 19 (a), 12 — 32 (b), 27 — 47 (c), 5 — 14 (d), and 7 — 19 (e). For comparison , we also show 
the crit i cal behavior of t h e resistivity in Si : P, a 3D system ( f ). Fig ures adapted from 
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We take the point of view that, at an empirical level, 
the experimental data presented in FigQ] represent a 
metal-insulator transition in 2D. It is, of course, possi- 
ble that the non-analytic behavior would ultimately be 
rounded out if it were possible to follow the physics to 
much lower temperatures than have yet been attained, 
but this would necessarily involve a crossover to new 
physics. However, the physics in the currently accessi- 
ble range of temperatures is important to understand in 
its own right. 

B. Strongly correlated highly metallic samples (r s 3> 1 and 

p < h/e 2 ): 

In this section we focus on samples with p <C h/e 2 . In 
this limit, there is a small parameter, 1/k-pt <C 1, which 
permits a well controlled perturbative expansion of phys- 
ical quantities. Because the predictions of Fermi liquid 
theory are sharp, discrepancies between experiment and 
theoretical expectations can be readily documented. 

1. Temperature dependence of p(T). 

As shown in Figs. Q] (a-e) and [21 in samples with large 
r s and p <C h/e 2 , the resistance p(T) is observed to 
increase with increasing temperature to a peak value 
p(T max ). Depending on which type of device and the 
value of n, the ratio p(T max ) / p{Tq) ranges from 2 to 10, 
where To ~ 25 mK is the lowest temperature at which 
resistivity measurements are readily carried out. 

This behavior has been observed in many 2D elec- 
tronic systems: Si MOSFETs, p- and n-GaAs quan- 
tum wells, p- and n-SiGe quantum wells, and in n-AlAs 
quantum wells. Even at the maximum, the resistance 
Pmax = p(T max ) is often smaller (and sometimes much 
smaller) than the quantum of resistance. For instance, 
the curves in Figf2] (b), (c), and (d) and the correspond- 
ing curves in FigQ] (b), (c), (d), and (e) are deep in the 
metallic regime. Moreover, generally it is found that 
fcflT max ~ Ef, where E-p is the bare Fermi energy. 

As a function of increasing temperature when T > 
T max , the resistance decreases, by as much as a factor 
of two (see FigJ2J), before ultimately starting to increase 



again at higher temperatures where, presumably, scatter- 
ing from thermally excited phonons starts to be impor- 
tant. Still deeper in the metallic regime {i.e. at smaller 
r s ), the resistivity is a monotonically increasing function 
of the temperature (see the lowest curves in FigQ] (a-e)). 

2. Parallel field magneto-resistance. 

As shown in Fig. O in Si MOSFETs, p-GaAs hetero- 
junctions, and SiGe quantum wells with large r s and p <C 
h/e 2 , the resistance p(T,B\\) for temperatures T < E-p 
exhibits a strongly positive magneto-resistance, increas- 
ing by as much as an order of magnitude as a function 
of .By, before saturating for Bu > B* ~ Ep/gps- Here 
Pb is the effective Bohr magneton, g is the gyromag- 
netic ratio, and B\\ is the magnetic field parallel to the 
film. In sufficiently thin samples, B\\ has little effect on 
the orbital motion of the electrons, so it can be viewed 
as coupling only to the electron spin, and therefore the 
magneto-resistance is directly a function of the degree of 
spin polarization of the electron liquid. 

In samples that are sufficiently close to the point of the 
metal-insulator transition, B\\ can even induce a metal 
insulator transition. This effect has been seen in Si- 
MOSFETs (see FiggJ). 

As shown in FigJS] (a) and (b) , the T dependence of p at 
low temperatures (T < T max ) is largely eliminated when 
the electron spins are polarized. Specifically, for magnetic 
fields B\\ > B*, the slope dp(T,B\\)/dT is reduced from 
its Bu = value - in some cases by as much as two 
orders of magnitude! This effect has been observed in Si 
MOSFET's and p-GaAs quantum wells. 

3. Magnetoresistance in a perpendicular magnetic field 

The experimentally observed behavior of 2D strongly 
correlated electron liquids in a perpendicular magnetic 
field B± can be quite complex, as shown in Fig[6l In 
part, this complexity reflects a combination of orbital 
and spin effects. 

At small B± metallic samples often exhibit relatively 
small negative magneto-resistance. At larger B±_ the 
magneto-resistance is typically large and positive, while 
at still larger B± it becomes negative again as the system 
enters the quantum Hall regime. 
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FIG. 2 Non-monotonic temperature dependence of the resistivity in a (100) Si MOSFET (a), p-GaAs quantum well (b), p- 
SiGe quantum well (c), and (111) Si MOSFET (d) deep in the metallic regime over an exten ded temperature range. The bare 
(no n-renormalized) F ermi tem peratures a re 7.5 K (a), 0.75 K (b), and 7 K (c). Adapted from llVlokashi and Kravchenkol |200g|) 
(a). bao et aU fcood ) (b), and lcolerideel U997J) (c). Panel (d) courtesy of R. N. McFarland and B. E. Kane. 



The interplay between the quantum Hall states and 
the behavior of the system at large r s and B± — ► has 
only been partially explored. The phase boundary in the 
n — B± plane between the quantum Hall and insulating 
regimes can be identified in experiment in at least two 
ways: a) Since the resistance is an increasing function of 
T in the quantum Hall phase (vanishing as T — > 0) and 
a decreasing function in the insulating phase, the phase 
boundary can be approximately identified as the points 
at which p xx is T independent, b) Since a xx vanishes as 
T — > in both the insulating and quantum Hall phases, 
but is non-zero in the critical regime, the phase boundary 



can be approximately identified as the points at which 
a xx has a local maximum at the lowest accessible T. At 
least when the field is not too small, these two methods 
produce essentially the same results, so the identification 
of the phase boundaries is unambiguous. 

As shown in the left-hand side panel of Fig. [Jj in a n- 
GaAs sample with relatively small r s , the phase boundary 
shifts up in energy as B — > 0, presumably reflecting the 
expected "float i ng" o f the extended states (jKhmelnitskii , 



1984 



Laughlinl . 119841 ) (i.e. for non-interacting electrons, 
if all the states at B± = are localized, then the lowest 
energy at which there is a localized state must diverge 
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FIG. 3 Giant magnetoresistance in a parallel magnetic field in a strongly metallic Si MOS FET (a), a 10 nm w i de p-GaAs 
quant u m well (b), a n-AlAs quan tum well (c), and a n -SiGe quantum well (d). Adapted from 
( 20od ). Ete"poortere eTal\ ( 20021 ) and Lai et ail 1 20051 ). respectively. 
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as B± — > 0). This "floating" can be seen even better in 
Fig. [8l where the energy of the lowest extended state in a 
strongly-disordered 2D hole system in a Ge/SiGe quan- 
tum well is observed to increas e by an order of magnitude 
as B±_ — > pilke et q.lhoool) . 

In contrast, in high mobility p-GaAs samples with 
large r s , as B± — > this phase boundary is observed to 
extrapolate to roughly the same value, n c , which marks 
the zero field metal-insulator transition in the same dc- 



of Fig. [7] The analogous phase boundary, which also 
extrapolates to n c , h as been traced in Si MOSFETs 



(jShashkin et all Il993l ). as shown in Fig. [9] While there 
are some notable differences in the way, within the quan- 
tum Hall regime, the various different integer quantum 
Hall phases terminate at low fields, in both p-GaAs and 
Si MOSFETs, the phase boundary between the quantum 
Hall and insulating phases clearly extrapolates to a finite 
zero field critical value. This further corroborates the 
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FIG. 4 Parallel magnetic field-induce d metal-insulator tran 
sition in a Si MOSFET. Adapted from lsimonian et all ( 19971 ) 



identification of critical point. 

At low fields, or elevated temperatures, and small r s , 
one naturally observes Shubnikov-de Haas oscillations of 
p{B±) in sufficiently clean 2D systems. Similar oscilla- 
tions have been observed in both the Si MOSFETs and p- 
GaAs quantum wells with large r s (see Fig. El). However , 



Gao et al 



(2003), 



for the large r s samples reported by 
these oscillations persist up to temperatures which are 
comparable with the bare Fermi energy. It is important 
to stress that FL theory not only predicts the existence of 
magnetic oscillations with period inversely proportional 
to the area enclosed by the Fermi surface, but it also 
predicts that the amplitude of these oscillations should 
decrease in proportion to exp[— 2irkBT/?iuj c \ where the 
cyclotron energy, u> c — eB/mc, in turn must be less than 



4. Spin magnetization of the electron gas at large r s 

The magnetization of the large r s metal is widely 
observed to exhibit a strong dependence on n, which 
reflects the increasing importance of electron corre- 
lations. For instance, the in-plane magnetic field 
of complete spin polarization, taken from magnetiza- 
tion and mag netocapacitance meas urements in clean 



Si-MOSFETs (Shashkin et al. 



20061 ). decreases signifi- 
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FIG. 5 The slope da/dT as a function of the parallel magnetic 
field in a Si MOSFET (a) a nd in a 10 nm wid e p- GaAs quan 
turn w ell (b) . Adapted from 
(2006), respectively. 
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n, as shown in Fig.[Tni Strong dep endences of the magne 



tization on n h a ve als o been seen (Vitkalov et al. (2001): 



Pudalov et al 



(2002): 



Zhu et 



((2003|); 



Vakili et al. 



(|2004l ); iLu et all ((2008)) in other types of devices with 
n near the critical density for the metal-insulator tran- 
sition. However, there are subtle, but important dif- 
ferences in the n dependences. For instance, in Si- 
MOSFETs, B* a ppears to extrapolate t o at a positive 
value of n = 



( Shashkin et all (|200l| ) 



pOOlh : IShashkin et all pOOtJ )). while in t he Si/ SiGe 



Vitkalov et al. 



cantly more strongly with decreasing n than docs Ep cx 



quantum well devices studied by ILu et all (j2008f ). the 
n dependence of B* over the accessible range of n was 
fitted by B* ~ n 13 , which vanishes only as n — > 0. 

Complex evolution of the linear susceptibility, x-i with 
n has also been observed, but in contrast with the non- 
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FIG. 6 Magnetoresistance of the strongly correlated metallic 2 D systems in a pe r pendi cular magnetic f ield in a Si MOSFET 
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(20(XJ) and 



Gao et al 



(2003), respectively. 
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Adapted from bultz et all jl99sl ). 



linear response (as parameterized by B*), the linear re- 
sponse can be strongly affected by rare regions of local- 
ized spins whi ch make a l a rge co ntribution to the low T 
suscep tibility (IPrus et all 120031) . Nevertheless, experi- 



ments (Shashkin et al 



20061 ) show that the n dependent 



trends of \ are roughly in agreement with the results 
extracted from measurements of B*(n). 



results from different systems, one vexing issue is to what 
extent the large values of p reflect strong effects of disor- 
der, as opposed to the intrinsic effects of correlations. We 
are primarily interested in data on the cleanest possible 
systems, where at the very least the correlation effects 
must have strongly enhanced the effects of the weak dis- 
order. 



C. Strongly correlated highly resistive samples, r s 3> 1 and 

p > h/e 2 

The properties of the 2DEG at large r s well on the in- 
sulating side of the metal insulator transition have been 
less completely explored, experimentally. In comparing 



Here, we focus on some astonishing experimental ob- 
servations in strongly correlated, highly resistive devices. 
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FIG. 8 A map of the extended states for a highly-disordered 
2D hole system in a Ge/SiGe quantum well. The open circles 
represent the positions of the extended state in the quantum 
Hall effect regime. The solid circles correspond to the position 
of the lowest extended state. Numbers show the value of 



pxyh/e 2 . Adapted from lHilke et all |2000|) 




B_ L (T) 

FIG. 9 Bands of the extended states (white areas) in a di- 
lute Si MOSFET. The shaded area corresponds to localized 
states. The i nset shows the exp e cted "floating" of the ex- 



tended states (Khmclnitskii 



from Shashkin et al 



1 199a) 
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Lauehlinl . Il984 :- 1. Adapted 
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FIG. 10 The critical magnetic field B* needed for complete 
spin polarization in a Si MOSFET as a functio n of the electron 
density. Adapted from Ishashkin et al. ( 20od l . 
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FIG. 11 Metallic (dp/dT > 0) temperature dependence of 
the resistivity in a 30 nm wide p-GaAs quantum well with 
p > h/e 2 (see four lower curves). 2D hole density is 6.1, 5.9, 
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1. "Metallic" T dependence in samples with resistivity 

p > h/e 2 

The conventional theory of localization in the strong 
disorder limit p ^> h/e 2 predicts that the electronic trans- 
port should be due to hopping conductivity, and hence p 
should diverge strongly (exponentially) as T — > 0. On the 



other hand, in p-GaAs samples with large r s , as shown 
in Figllll the low temperature resistivity for a range of 
n in which p is up to 3 times larger than h/e 2 some- 
times exhibits a "metallic" temperature dependence, i.e. 
p increases with increasing temperature. Similar behav- 
ior has also been observed in Si MOSFETs (see the two 
middle curves in Fig[T] (a)). 
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FIG. 12 Parallel field magnetoresistance in a Si MOSFET at 
different electron d ensities across the m etal-insulator transi- 
tion. Adapted from Mertes et al. ( 19991 ). 



2. Magneto-resistance for r s 2> 1 and p > h/e 2 . 

The magneto-resistance in a parallel magnetic field Bu 
on the insulating side of the metal-insulator transition is 
still large and positive, even in samples with p as high 
as 2 Mil/square. This can be seen in the data from Si 
MOSFET devices shown in Fig. [12] At large enough B\\ 
the magnetoresistance saturates in a way similar to that 
in metallic samples. 

In a perpendicular magnetic field, B±, samples with 
p(B = 0) 3> h/e 2 exhibit rather complex magnetoresis- 
tance. As can be seen in Figs. [I3](a-c), the resistance first 
increases strongly for small B±, sometimes by orders of 
magnitude. At higher B±, the magnetoresistance turns 
negative and eventually the system exhibits giant resis- 
tance oscillations with deep minima at magnetic fields 
corresponding to Landau level filling factors v = 1 and 2. 
The sam e general fe atures have been observed in p-GaAs 
samples (jGaol 120031 ). as shown in Fig[l3] (d)). 

That the positive contribution to the magneto- 
resistance is largely a spin effect can be seen in Fig. [T4l 
where, in the presence of a strong in-plane field, Bu ~ B* , 
the resistance at B± = is much larger than for Bu = 0, 
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FIG. 13 Magnetoresistance of a highly resistive dilute 2D gas 
in a perpendicular magnetic field in a Si MOSFET (a-c) and 
in a 30nm w i de p- GaA s quantum well (d). Adapted from 
D'lorio et all (jl99uh and baol (J200J), respectively. 



but then, as a function of increasing B±, the magneto- 
resistance is everywhere strongly negative all the way to 
the quantum Hall regime. 



D. Drag experiments on double-layers with r s 3> 1 and 

p < h/e 2 

Additional information concerning correlation effects 
can be obtained from measurements of the "drag" resis- 
tance in a system of two 2DEG layers which are electri- 
cally unconnected. Current / is passed through the lower 
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FIG. 14 R XX {B ± ) in a silicon MOSFET in the presence of 
a parallel field Bu — 3.4 tesla used to suppress the metallic 
behavior (solid symbols). For comparison, the magnetoresis- 
tance in the absence of a pa rallel magnetic fiel d is sh own by 
the solid line. Adapted from iKravchenko et al\ (| 19981 ). 



(active) layer and the voltage Vjj is measured on the up- 
per (passive) layer. The drag resistance is defined to be 
the ratio pu = Vb/J. 

Fo r relati vely small r s , measurements (jGramila et al 



1993 



19941 ) of p D (T) in double layer 2DEGs are in 



quali t ative agreement with F e rmi l iquid theory (jPrice 



1983 



Zheng and MacDonaldl . I1993T) . Specifically, the 
drag resistance is small, proportional to (T/E-p) 2 and 
to (kF<i)~ ad , where fcp is the Fermi momentum, d is the 
spacing between the two layers, and typically ad = 2 or 
4, depending on the ratio of d/i. 



However, experiments (jPillarisettv et all 120031 ) on p- 
GaAs double layers with r s ~ 20 — 30 yield results which 
differ significantly from what would be expected on the 
basis of a simple extrapolation of the small r s results. 
These experiments were performed on samples with small 
resistances p ~ 0.05 — 0.1 h/e 2 , in which quantum inter- 
ference corrections to the Drude conductivity are pre- 
sumably insignificant. The experimental data shown in 
Fig. [15] reveal the following features: 

• The drag resistance in these samples is 1-2 orders 
of magnitude larger than expected on the basis of 
a simple extrapolation of the small r s results. 

• Whereas in a Fermi liquid, pd{T) ~ T 2 . in large r s 
devices po ~ (T) aT where the temperature expo- 
nent exhibits non- Fermi liquid values 2 < ctT < 3. 
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FIG. 15 In- plane magnetotransport data for p m = 2.15 x 10 10 
cm" 2 at T = 80, 175, 250, and 400 mK. (a) Inset: p vs By . 
Main Plot: Data from inset normalized by its B« = value, 
(b) Inset: Corresponding data for po vs Bu. Main Plot: 



Data from inset normalized by its B\\ 
from lPillarisettv et al\ ([200J). 



II- 

value. Adapted 



• At low temperature, pd{T, -By) increases as a func- 
tion of Bn by a factor of 10-20 and saturates when 
£>H > B* . A parallel magnetic field also appears 
to strongly suppress the temperature dependence 
of po{T). In the presence of a non-zero By, the 
value of olt{B\\) decreases with Bu and saturates 
for B\\ > B* at a value which is significantly smaller 
than the FL value ut = 2. 

• The T and especially the Bn dependences of 
Pd{T, B\\) and the resistivities of the individual lay- 
ers p(B\\ , T) look qualitatively similar to one an- 
other, which suggests that both have a common ori- 
gin (see Fig. I15[) . (In a Fermi liquid, p(T, Bu) is pri- 
marily determined by the electron-impurity scat- 
tering, and pd{T, B\\) by the inter-layer electron- 
electron scattering, so there is no a priori reason 
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for their T and B\\ dependences to be similar.) 

E. Comparison with small r s devices 

The anomalies discussed above have been observed in 
samples with large r s . In relatively smaller r s high mobil- 
ity devices ( i.e., r s ~ 1), the observed behavior is much 
more in line with the expectations of FL theory, modi- 
fied by weak interference effects. For example, there is 
a smooth finite temperature crossover as a function of 
disorder from a "weak localization" regime for p < h/e 2 
to variable range hopping behavior for p > h/e 2 ; p is 
weakly (logarithmically) dependent on T for p < h/e 2 
and T < Tj?; the parallel field magneto-resistance is 
weak (logarithmic) and positive for p < h/e 2 ; the re- 
sistance diverges strongly (exponentially) with decreas- 
ing T whenever p > h/e 2 ; the d rag resistan c e is s mall 



Ando et al 



( 1982ft for 



compared to p (see, for example, 
a review). Moreover, in the presence of a perpendicu- 
lar magnetic field, small r s devices exhibit "levitation of 
the delocalized states", i.e., the phase boundary between 
the quantum Hall and the insulating phases tends toward 
ever higher values of the density as B± — > 0, as shown 
on the left hand side panel of Fig. [7] and in Fig. [H 

II. THEORETICAL CONSIDERATIONS 

A. Good 2D "metals" with r s < 1 and k f £ > 1 

The properties of interacting electronic systems cru- 
cially depend on the dimensionless strength of the in- 
teractions, r s , and the strength of the quenched disor- 
der, which is parameterized by the dimensionless quan- 
tity kp£, where k-p is the electron Fermi momentum and 
t is the electron elastic mean free path. The theory of 
pure (kp£ ^> 1) and weakly interacting (r s <C 1) elec- 
tron liquids is under good theoretical control. In this 
case, Boltzmann theory yields a good first approxima- 
tion to the transport properties, and any Fermi liquid 
corrections to the bare electron mass and to the den- 
sity of states are small. Moreover, the screening radius 
A sc = Qjg/4, is much larger than the average distance 
between electrons: X sc n 1 / 2 = (4y / 7r) _1 r^ 1 » 1. 

In this limit, to a first approximation, the transport 
scattering rate can be calculated by perturbation the- 



ory. Specifically p is proportional to the rate of trans- 
fer of momentum from the electrons to the lattice. At 
high temperatures the resistance of the 2D electron gas 
is determined by the clcctron-phonon scattering. At tem- 
peratures low compared to the effective Dcbye temper- 
ature, T p h ~ hck-p, where c is the speed of sound, the 
electron phonon scattering is no longer significant. (In 
some cases, a crossover temperature, T p h, below which 
electron-phonon scattering is unimportant, can be read- 
ily identified; for example, in Figs. [2j T p h > Ep is roughly 
the point at which p(T) has a minimum.) 

Since in all cases of interest here, the Fermi momen- 
tum is much smaller than the reciprocal lattice vec- 
tor, electron-electron scattering conserves the total quasi- 
momentum to very high order, and therefore does not 
contribute directly to the resistance. (This is different 
from metals with large Fermi momenta where electron- 
electron Umklapp processe s determine the T -dependence 
of the resistance at low T (jAbrikosovl . 119881 ).) Therefore 
the low temperature resistance of the system is due to 
electron-impurity scattering. 

Consider, for example, the case in which the dominant 
low temperature scattering is from the Coulomb poten- 
tial due to a random distribution of charged impurities. 
Naturally, the potential of an impurity is screened by the 
2DEG, itself. Two generic aspects of the screening are 
that it causes the potential to fall more rapidly when the 
distance from the impurity exceeds the screening length 
and it induces a rapidly oscillating component of the 
screened potential corresponding to the Fricdel oscilla- 
tions in the density. As is well known, scattering by an 
unscreened Coulomb potential produces singular forward 
scattering, leading to an infinite cross section. However, 
the transport scattering rate weights large angle scatter- 
ing more heavily In 3D, the transport scattering cross 
section is still logarithmically divergent, so even when 
the screening length is long, it plays an essential role in 
cutting off this divergence. In 2D, however, the transport 
scattering rate is finite, even in the limit \ sc — > oo, where 
in the Born approximation it is given by 

h _ (27r) 2 e 4 
~ ~ Ep 



-N, 



l {i + o([k F x sc }- 2 )} 



(1) 



where is the concentration of impurities. Moreover, to 
the extent that quantum inference between different scat- 
tering events can be neglected (Boltzmann transport), 
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the conductivity is related to t, the Fermi velocity, vp, 
and the density of states, v, according to the Drude for- 
mula 



In infinite samples L in Eq. [4] should be reinterpreted as 



e l v F rv/2. 



(2) 



Here, p approaches a constant value as T — ► 0. A parallel 
magnetic field, which we assume acts only on electron 
spins, changes the spin degeneracy at the Fermi level. In 
the limiting cases, Bu = 0, there is a spin degeneracy 
g c = 2, while for Bu 3> B* , g c = 1. Since E-p oc vp oc 
g~ , and v oc g c , we conclude that for r s « 1, 

P « .9c- (3) 

This means that the resistance decreases by a factor of 
2 as .By increases from to B*. By the same token the 
scale for the T dependence of p(T) is set by Tp = Ep/ks] 
since at temperatures large compared to Tp, t oc v f oc T, 
it follows that p decreases as T increases in proportion 
to T~ 2 . The only change in the T dependence of p in 
the presence of B\\ > B* is produced by the factor of 2 
increase in Tp. Similar considerations (which we will not 
review explicitly) lead to the conclusion that the drag re- 
sistance in double layers is small in proportion to T 2 / Ep 
and consequently po decreases as B\\ increases. 

This simple theoretical description in the limit r s «l 
and kpl 3> 1 gives a good account of the transport in 
clean, high density devices, but is in drastic qualitative 
disagreement with all the experimental results on the 
transport properties of the large r s > 1 devices presented 
in Section 1. 



B. Weak localization corrections and the theory of 2D 
localization 

In 2D, the interference between multiple scatter- 
ing processes, i.e., the "weak localization correc- 
tions" to the conductivity , dive r ge at T — > 0, an d 



oo (Abrahams et al 



1979 



Gorkov et al 



1979). 



Here L is the samp l e size . (For a review, see 
Lee and Ramakrishnanl (|1985l ).) This divergence, how- 
ever, is only logarithmic, and therefore for kpl 3> 1, 
and at accessible temperatures, these corrections are 
small in comparison to the zeroth order conductivity 
G = (h/e^p- 1 =G + SG, where G ~ kpl and 



the phase breaking length = 'Dtq, where D = vp£ 
is the diffusion constant and ~ T p is the phase- 

breaking rate, in which p depends on details of inelastic 
scattering processes. 

C. Interaction corrections 

Impurities in a metal create Friedel oscillations of the 
electron density. Due to the electron-electron interac- 
tions, the quasiparticles in the metal are scattered not 
only from the impurity but also from the modulations 
of the electron density. The interference between these 
two processes gives rise to corrections to the Drude resis- 
tance. These corrections are interesting because they are 
non-analytic functions of T and B\\ , so at small enough 
temperatures they dominate the T and Bu dependences 
of the resistance. 



In the diffusive limit (L T = y/D/T > I) the 



cally divergent as Lt - 


-> oo {£> 


dtshuler et al. (1980): 


Finkelshtein (1983): 


Fin] 


celshtcm 


(1984ah: Finkelshtein 


(1984bl): Castellani et al. 


(1984): 


Altshuler and Aronov 


(1985): FinkelsteirJ ( 


199(1)): 



5G = - 



2ir 2 { 



1 + 3 



1 - 



ln(l + Fp) 

Fn 



}\n(L T /£), (5) 



where Fa < is an interaction constant in the triplet 
channel. Note that this same interaction parameter is 
also responsible for an enhancement of the spin suscep- 
tibility, x — Xo/0- + Fa). For r s <C 1 (in which limit, 
|i*o| -C 1), Eq. Ogives (up to a numerical factor) the same 
(negative) correction to the conductivity as Eq. [4] Note, 
however, that this correction has the opposite ("metal- 
lic") sign when ln(l + F )/F > 4/3. 

At higher temperatures, when Lt "C £, the leading 
interaction effect involves the interference between a sin- 
gle electron scattering from an impurity, and from the 
Friedel oscillations in the neighborhood of the same im- 
purity. The interference corrections to the Drude for- 
m ula in this so called "ballistic regime" were considered 

by 



Zala et al 



(2001 



5G 



ln[L/4 



(4) 



The result is that in an intermediate interval of tem- 
peratures, E-p ^> T ^> h/r, there is a T-dependcnt cor- 
rection to the conductance which is is linear in T. It 
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consis ts of the sum of " singlet" and "triplet" contribu- 
2001allbh : 



tions ( Zala et 



G(T)-G(0)ocG -f- 
&p 



1 



Fa 



l + F 



(0) 



The factor of 3 in the second (triplet) part of Eq.[6]refiects 
the existence of 3 channels in the triplet part on electron- 
electron interaction. (Eq. is written for the case in 
which there is no valley degeneracy, as in GaAs, but it 
is readily generalized to the valley degenerate case. For 
example, in the case of Si MOSFETs, with a two-fold 
valley degeneracy, this factor becomes 15.) For repulsive 
interactions Fq is negative with a magnitude which, for 
small r s , is proportional to r s . Thus, at small r s , p(T) 
decreases linearly as T increases. 

Of course, most experimental realizations of the 2DEG 
have r s > 1 . At a phenomen ological level, it is possible to 

extrapolating to stronger 



imagine (jZala et 



2001 



33) 



interactions in the spirit of Fermi liquid theory, in which 
case for Fq < —1/4, a linearly increasing resistance would 
result. The range of validity of this sort of extrapolation, 
which is clearly sensible to some degree, is a unresolved 
theoretical issue in the field, to which we will return. 



D. The theory at r s S> 1 

In an ideal MOSFET, where there is a metallic ground 
plane displaced by a distance d from the 2D electron 
gas, the electron potential energy changes its form as a 
function of d. For r -C d electrons interact via Coulomb 
interaction and V(r) ~ e 2 /r. For r 3> d one has to take 
into account the interaction between electrons and their 
images in the ground plane and hence V(r) ~ e 2 d 2 /r 3 . 

The mean-field phase diagram for this problem is 
shown schematically in Fig ll6l The considerations that 
produce the re-entrant character of transition line as a 
function of n are relatively simple. The kinetic energy 
of the 2D electron liquid scale with the electron density 
as ~ n/m. For Coulomb interactions, the typical inter- 
action strength scales as ~ e 2 ™ 1 / 2 , so the ratio of the 
potential to kinetic energy, r s oc n -1 / 2 , decreases with 
increasing density. Therefore, in the range of densities 
n ^> dr 2 , increasing density always favors the FL. Con- 
versely, for n <C d~ 2 , the interactions between electrons 
are dipolar, and hence effectively short-ranged, i.e. the 



typical interact scales as ~ e 2 d 2 r?l 2 . Thus, in this range 
of density, r s oc n 1 / 2 decreases with decreasing density, 
and hence decreasing density favors the FL, reflecting 
the same trends as does 3 He. One implication of this 
analysis is that there exists a critical value of the ratio, 
d/ap = ot c , such that for d < a c ap, there is no WC phase 
at any density. In this case there is a largest achievable 
value of Max[r s ] ~ d/a* B . 

Fixed node Monte Carlo simulations 



(jTanatar and Ceperlevl . 119891 ) performed under the 
assumption that there is a direct FL-WC transition yield 

(c) 

the large critical value = 38 ^ 1. By analogy, we 
think that it is likely that a c ~ is also large compared 
to 1, although as far as we know, this issue has not been 
addressed using any quantitative methods. 

The very existence of a "highly correlated" 2D elec- 

(c) 

tron liquid in the range 1 < r s < derives from the 

(c) 

large value of r$ . From this point of view, it is pos- 

(c) 

sible to consider the case r s < r$ , and still treat r s 
as a large parameter. If we do this, we see that the 
highly correlated 2D electron liquid has (at least) three 
characteristic energy scales: 1) E-p, or more properly 



( Andreey 



19781 Il979i ). the rcnormalizcd Fermi energy, 



E F < E-p, which contains a rcnormalizcd mass, 2) the in- 
teraction energy V = r s Ep, and 3) the plasma frequency 
ftp ~ \[EpV = ^frlEp. For large r s , these energies are 
quite distinct, with V > ftp > Ep. As a consequence 
of the existence of this hierarchy of energy scales there 
arc 4 distinct temperature intervals where the electron 
liquid behaves differently. We will discuss these inter- 
vals below: T < Ep where the system is in the Fermi 
liquid state, Ep < T < rip where the system is a non- 
degenerate strongly correlated and still highly quantum 
liquid, flp < T < V where the system is a highly corre- 
lated classical liquid, and T > V where the system is a 
classical electron gas. 



E. The Fermi fluid to Wigner crystal transition in the 
absence of disorder 

The classic studies of the Fermi fluid to WC tran- 
sition were carried out under the assumption that, 
at T = 0, there is a direct first order transi- 
tion where the ground stat e energies cross. How- 



ever, it has been shown (jLorenzana et all (|2002j); 
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Microemulsion phases. 



FIG. 16 Schematic phase diagram for 2D electrons in a MOSFET with no disorder. 



Spivakl (120031); ISpivak and Kiyelsonl (12004); Ijamei et al 



(2005); ISpivak and Kivelsonl |2006l) 



Giuliani et al 



20071 ); iBiskup et all (|2007|)) that for the case of a rel- 



atively long range potential 



V 



1 > x > 3, 



(7) 



a direct first order liquid-crystal transition is forbidden in 
2D. Instead, either the freezing transition is continuous 
(whic h is unlikely, for well known reasons (jBrazovskii . 
19751 )). or between these two phases there must oc- 
cur other phases, which we have called "microemulsion 
phases" . 

The largest theoretical uncertainty, here, is that no re- 
liable estimates exist concerning the width in density of 
these novel phases. The part of the phase diagram where 
these phases should exist is indicated qualitatively by the 
shaded region in Fig |16l A rich variety of phases, includ- 
ing bubble and stripe phases, are expected in this region. 
The existence of a stripe phase, for example, would be 
detectable macroscopically through a spontaneously gen- 
erated resistance anisotropy. 

At the moment there is no direct experimental evidence 
of the existence of such phases. This is either due to the 
fact that, for some reason, the range of densities in which 
they exist is extremely small, or because currently avail- 
able samples are not pure enough to reflect the physics of 
the zero disorder limit. In any case, the existence of such 
phases, which are very different from both the FL and 
the WC, significantly complicates the theory of strongly 
correlated disordered electronic systems. 



F. Theoretical considerations concerning the 
metal-insulator transition in 2D with disorder 

1. Anderson localization 



It has been suggested ([Abrahams et all Il979l ) that in 



absence of electron-electron interactions the logarithmic 
derivative 



dhxG 



/3(G) 



(8) 



rilnL 

is a function of G alone. According to Eq. [H (3 ~ —1/G 
as G — > co. The fundamental result of the single particle 
theory of localization follows from this: even weak dis- 
order leads to insulating behavior at T = 0, and hence 
there can be no metal-insulator transition. On the face 
of it, this statement is inconsistent with the experiments 
summarized in the previous section. 

2. Transition in the presence of spin orbit scattering 

Since the existence of a metal insulator transition is 
an important issue, it is worth considering the single 
particle problem in the presence of spin orbit coupling, 
where the theory predicts a zero temperature transi- 
tion between an " ideal metal" ( G = o o) and an in - 
sulator (G = 01 fflikami et all |l980); Undol |l989T ); 



Evangelou and Zimanl (|1987[ )1. Analytically, this follows 
again from a pcrturbative analysis of the interference cor- 
rections: 8G = +j\n(L/£ so ) for L large compared to 
the spin-orbit scattering length, L £ so > £. In other 
words, so long as G at length scale ~ i so is large com- 
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pared to 1, at low enough temperatures that > £ so , 
the conductivity is an increasing function of decreasing 
T - "antilocalization." Under the same assumption of 
one-parameter scaling as in Eq. [5J this implies that 
ft + 1/G as G — > oo, which in turn implies the sta- 
bility of a perfectly metallic phase with p = 0. On the 
other hand, the stability of an insulating phase for strong 
enough disorder is not debatable, i.e. ft < for small 
enough G. Moreover, the validity of the one-parameter 
scaling in this problem and the existence of a metal in- 
sulator transition have been tested in numerical studies 



(Asada et al 



2004) 



While this is the simplest model system in which a 
2D metal insulator transition exists, it is unlikely that 
it is relevant to the experiments discussed in Section 1. 
Specifically, the spin orbit interaction in n-silicon MOS- 
FETs is quite weak, and hence unlikely to be important 
in the accessible range of temperatures. 



3. Scaling theories of a metal-insulator transition with strong 
interactions 

Even in the absence of spin-orbit coupling, assuming 
the theory remains renormalizable, the ft function is fun- 
damentally modified by strong interactions. A key idea 



transition ( 


McMillan 


( 


1981 


); 


finkelstein 


( 


1990 


): 


Dobrosavlievic et al. 




1997 


); 


Punnoose and Finkelstein 


(2001 


); 


Punnoose and Finkelstein 




(200fi|)) is 


that the 



logarithmically diverging interaction correction to the 
conductivity in Eq. [5] can have the opposite sign 
and larger magnitude than the single particle weak 
localization corrections in Eq. |4] when extrapolated 
to large values of r s where \Fq\ is no longer small. 
Moreover, even where perturbation theory is still valid, 
i.e., where G 1 and r s <l. the value of Fq is subject 
to renormalization, which in the case of weak Coulomb 



interactions, \Fq\ <C 1, is of the form (jFinkelsteinl . I1990T ) 
dF 



dln(L/l) 



1 



(9) 



Since Fq is negative, this means that its magnitude grows 
with increasing L. It is therefore possible to imagine (in 
the spirit of the renormalization group, RG) that even 
where the bare interactions are weak, at sufficiently long 



scales, Fo grows until \Fq\ < 1, at which point the ft 
function would change sign. 

If a metal insulator transition occurs, it must occur 
at finite (dimcnsionless) interaction strength, F* , and 
finite conductance, G* . (Here, F may denote the 
strength of one or several interactions.) A priori, the 
properties of such a fixed point cannot be computed 
perturbatively. Either, the existence and character of 
such a fixed point must be inferred by extrapolating per- 



streneth (Castellani et al 




1984; 


Chamon and Mucciolo. 


2000; Finkelstein!. 199ol Navak and Yans 




2003: 


Punnoose and Finkelstein!. 


2001 


. 2003). or it must 


simply be conjectured ( 


Dobrosavlievic et al.. 


199' 


? ) on 



phenomenological grounds. One problem with this 
scenario is that for G ^> 1 , if the RG procedure is taken 
literally, the system evolves to a low temperature state 
where Fo — * — 1, which in addition to changing the sign 
of the ft function, implies the existence of a magnetic 
instability, which surely must affect the physics. More 
generally, there exists no clear qualitative picture of 
what happens to the system when the parameter Fq is 
significantly rcnormalized. 

There are a number of experimentally relevant conse- 
quences of this scenario. Most importantly, it implies the 
existence of a true quantum phase transition between a 
(possibly ideal) metal and an insulator, with all the im- 
plied quantum critical phenomena. Secondly, valley de- 
generacy figures as an important factor for the character, 
and po ssibly even the existence of a metal-insu lator tran- 



sition (Punnoose and Finkelstein 



2001 



20051) . Thirdly, 



the theory predicts the existence of a peak in the tem- 
perature dependence of p(T). For systems with kp£ ^> 1, 



perature ( 


Finkelstein 


2001, 


2005h 



1990: 



Punnoose and Finkelstein 



Tmax ~ (vf/#) exp(-aG) 



(10) 



corresponding to the length scale, £(T max ), at which ft 
changes sign. Here a is a constant of order one, which 
depends on the strength of the bare interactions. For 
T < T max , the system exhibits "metallic" T dependence. 
The theory also predicts an increase of the resistance at 
small values of B\\ . However when the conductance of the 
system is large, G>1, both these effects become very 
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small and manifest themselves only at exponentially low 
temperatures (see Eq. QI 



G. 2DEG in the presence of a perpendicular magnetic field 

In considering the phase diagram of the 2DEG in the n 
- B± plane with a fixed strength o f disorder, ther e are a 
number of asymptotic statements (jKivelson et all Il992n 
that can be made with confidence on theoretical grounds: 

• Firstly, for low enough density n < n m i n and for 
any strength of magnetic field the system must be 
insulating; in the absence of disorder, this is a con- 
sequence of Wigncr crystallization, and in the pres- 
ence of disorder, it can be traced to the strong ten- 
dency of dilute electrons to be strongly localized by 
disorder. 

• For any fixed density, and high enough magnetic 
field, n(f>o/B± <C 1, it is similarly straightforward 
to show that the system must always be insulating. 
(Here, cf>o = hc/e is the quantum of flux.) 

• Most interestingly, in the limit that both n and B± 
are large, but with n<j)o/B± ~ 1, there necessarily 
exist robust quantum Hall phases. To see this, note 
that in this limit, Tilo c is large both compared to the 
strength of the disorder and the electron-electron 
interactions, so inter-Landau level scattering can be 
treated perturbatively. (Here, u) c = eB±/mc.) In 
this limit, the existence and stability of integer and 
fractional quantum Hall states is well established. 

For present purposes, we will neglect the interesting 
complexity associated with the various distinct quantum 
Hall phases, and simply discuss the considerations that 
determine the shape and topology of the curve, n*(B±), 
that encloses the regions of quantum Hall phases in the 
phase diagram (see Figs (8] and [9|). From the three obser- 
vations above, it follows that 



n*(B ± ) ~ BJ(<t> ) as B ± 



(11) 



In the absence of electron-electron interactions and spin- 
orbit scattering, we know that all states at B± = are 
localized. This requires that n*(B±) — > oo as B± — > 0. 
Within the single-particle theory, this constraint was ac- 
counted for by the notion of "le v itation" of the d elocal- 
ized states (IKhmelnitskiil . Il984t iLaughlinl . |l984). The 



idea here is that in the large field limit, when lo c 1/r, 
there is a single delocalized energy level at the cen- 
ter of each Landau level, but that when at small field, 
u> c ~ 1/r, the delocalized states initially associated with 
each Landau level remain distinct, but begin to "levi- 
tate" to higher energies. Thus, as B± — > 0, the energies 
of the delo calized states dive rge. More recent numeri- 
cal studies (jSheng et aZ.I . l200lh have suggested that, even 
for non- interacting electrons, the fate of the delocalized 
states may be a more complex issue. However, clearly, so 
long as there is only a single, insulating phase at B± = 0, 
the divergence of n*(B±) as B±_ — > is inescapable. 

Conversely, if n*(B±) tends to a finite value, n c = 
n*(0), as B± — > 0, it probably implies the existence of 
a zero field phase transition at this value of the density. 
Since n* is unambiguously a critical density marking the 
point of a quantum phase transition, tracking it to the 
zero field limit is a very promising strategy for distin- 
guishing a crossover from a true phase transition. 

Comparing these theoretical expectations with the al- 
ready discussed experiments, we see that in the large r s 
p-GaAs and Si MOSFETs (the right hand side panel in 
Fig|7]and Figj9]) , n*{B±_) clearly extrapolates to a finite 
value as B± — > 0. Conversely, in the small r s n-GaAs 
device (left hand side panel in FigJ7]and Fig[5]), clear ev- 
idence of levitation is seen, in the sense that n*(B±) has 
a pronounced minimum at a value of B± which arguably 
corresponds to uj c t ~ 1, and then grows strongly as B± 
is further reduced. 



H. Microemulsions of WC and FL in the presence of weak 
disorder 

I. Effects of weak disorder at T = 

Starting from the clean limit, treating the interactions 
as strong and the disorder as a small perturbation leads 
to a rather different perspective on the metal insulator 
transition. It is important to recognize, however, that 
even weak disorder is always a relevant perturbation in 
2D. If we ignore the existence of microcmulsion phases in 
the zero disorder limit, then in 2D disor der always rounds 



the fir st order liquid to WC t r ansiti on ( Imry and Worti 



(|1979l ); lAizenman and Wehrj (| 19891) : 



Baker 



(11993|)), re- 



sulting in a phase in which islands and continents of 
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pinned WC coexist with seas and rivers of Fermi fluid. 
If, instead, we start with a spatially ordered pattern of 
alternating WC and Fermi fluid, characteristic of a mi- 
croemulsion phase, then disorder destroys the long-range 
order, again resulting in a phase which consists of a dis- 
ordered mixture of WC and Fermi fluid regions. Though 
the two scenarios differ greatly for weak disorder in the 
degree of local organization of the coexisting regions, at 
long distances they are difficult to distinguish. 

Since at T = 0, a WC does not conduct, this picture 
calls to mind a percolation type metal-insulator transi- 
tion. Of course, this picture does not take into account a 
variety of possibly important effects, including interfer- 
ence effects which could ultimately turn the percolativc 
metal into an insulator. However, in the limit in which 
the characteristic sizes of the islands and continents are 
large, classical percolation could be a good description 
over a significant range of T. 

2. Low T thermal physics: the Pomeranchuk effect 

The Fermi-liquid to crystal transition takes place at 
r s 1. In the Fermi liquid state, the entropy density, 

Si 



crgy between spins in a Wigncr crystal is exponen- 
tially small, J ~ exp[— Q ^/ rJl where a is a numbe r 



FL 



n 



is small for temperatures T <C Ep. In con- 
trast, for T large compared to J, the spin exchange in- 
teraction between localized particles in the crystal state, 
the spin entropy of the crystal Sc ~ nln2 is relatively 
large and temperature independent. Since for large r s , 
J is e x ponen tially smaller than Ep (See, for example, 
Rogerl (|1984[ ).). there exists a broad range of tempera- 
tures, Ep > T > J, in which the crystal is the high- 
entropy state, > Sfl- In this range of temperatures 
and near the point of the crystal-liquid transition the sys- 
tem tends to freeze upon heating! Of course, at ultra-low 
temperatures, T <C J, or at high temperatures, T > Ep 
the entropy of the crystal is smaller than that of the liq- 
uid, so the Pomeranchuk effect disappears. 

The tendency of Fermi systems to freeze upon 
heating, known as the "Pomeranchuk effect," was 
originally dis covered in the framework of t h e the - 
ory of 3 He (jPomeranchuk 



1950; 



Richardson 



1997). 



where it has been confirmed experimentally. (See 
the phase diagram of 3 Hc in Fig fT7| In this regard, 
the Wigncr crystal of electrons and the crystal of 
3 He atoms are similar. Specifically, the exchange en- 



of order one (IRogeri (119841): IChakravartv et all (119991 ); 



Voelker and Chakravartvl (|200lh V Thus at low tempera 



tures the evolution of the random microemulsion phases 
as a function of T and B\\ is dominated by the en- 
tropy changes associated with the sp i n degrees of freedom 
ISpivakl . l2003t Ispivak and Kivelsonl . boO^booeh . 

In the case of the two phase coexistence at low temper- 
atures, the fraction of the system which is locally Wigncr 
crystalline, /wd increases linearly with T. This behav- 
ior will dominate the T dependence of many important 
physical properties of the system at low T. 

Similar considerations govern the Bm dependence of the 
phase diagram. Since the spin susceptibility xwc 3> Xfl, 
the corresponding magnetization il%c 3> Mfl at small 
Bu . Since the free energy of the system contains the term 
—MBh , there is a £?||-induccd increase of the WC fraction 
over a wide range of circumstances. 

Since the resistance of the WC and FL are very differ- 
ent, the Pomeranchuk effect is one of the most directly 
testable features of the microemulsion phases. Where the 
FL is the majority phase, it leads to the robust predic- 
tion that p(T, £?||) increases with T and B\\ at low T, 
reflecting the purely thermodynamic fact that /wc in- 
creases as T and B\\ increase. However, for Bu > B* , 
the T and Bu dependences of p{B\\) are quenched since 
there is no spin-entropy in the fully polarized system, so 
fwc(T,B\\) no longer depends strongly on these quan- 
tities. This effect can, in principle, produce arbitrarily 
large fractional changes in the low T resistivity. Even on 
the insulating side of the transition, where at T = a 
majority of the system is Wigner crystalline, a tempera- 
ture induced increase in /wc can result in an increase of 
p with increasing T, in a limited range of temperatures, 
even though p > h/e 2 . 

Many other, more detailed aspects of the T and 
Sii dependences of p depend on the type of disor- 
der and other "details." While some progress has 
been made in understanding the dynamical properties 
of microemulsion phases i n the presence of weak disor- 
der (jSpivak and Kivelsonl . 120061 ). no satisfactory and/or 
quantitative theory of disordered microemulsions is cur- 
rently available. 

Of course, at high enough temperatures, where the 
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electron liquid is non-degenerate, the WC always is the 
lower entropy phase. The Pomeranchuk effect reflects, 
more than anything, the low entropy of the degenerate 
Fermi liquid at T < Ep . The melting temperature of the 
WC is determined by the (possibly strongly renormal- 
ized) Fermi energy of the "competing" fluid. Above this 
temperature, the microemulsions give way to a uniform 
non-dcgcncratc electron fluid. 

3. Crossovers at higher temperatures (T > Tp) and large r s 

In the absence of disorder, Umklapp scattering, and 
electron-phonon scattering, the long-wave-length proper- 
tics of the electron fluid are governed by hydrodynamics, 
rather than by the Boltzmann equation. In this limit, 
we need to talk about the viscosity of the fluid, 77, rather 
than the conductivity, which is in any case infinite. For 
gentle enough disorder, we can think of the disorder as 
defining some form of effective medium through which 
the otherwise hydrodynamic fluid flows. This is a trans- 
port regime, which is not often recognized, in which the 
hydrodynamic healing length (e.g., the electron-electron 
mean-free path, £ c _ c ) is short compared to the distance 
between impurities, or the length scale over which the 
disorder potential varies. In this case, 



p OC X] 



(12) 



where the proportionality constant is in general compli- 
cated function of the strength and character of the disor- 
der potential, and r\ is the viscosity of the electron fluid 



in the absence of disorder. The conditions for the validity 
of this equation are generally violated at T — > 0, where 
£ c - c — > 00, but for sufficiently strong interactions, it can 
be satisfied down to moderately low T. 

Here, we make a few remarks concerning the viscosity 
of the uniform fluid at T > Tp. As mentioned previ- 
ously, for large r s , there is a broad range of temperatures 
in which Tm ~ If < T < V, where the electron fluid, 
although non-degenerate, is still highly correlated. (Here 
Tm is the melting temperature of the Wigner crystal.) 
It is generally observed that the viscosity of highly cor- 
related fluids is a decreasing function of increasing tem- 
perature. This observation, combined with Eq. [12] and 
the low temperature increase in resistivity produced by 
the Pomeranchuk effect, implies that for metallic samples 
with large r s , there should generally be a maximum in 
the resistivity at T ~ Tp. 

Looking at this problem more carefully, there arc, as 
noted in Section Hl.DI two distinct ranges of T to be con- 
sidered. Beyond this, our theoretical understanding of 
the viscosity of correlated fluids is crude at best. We thus 
rely on the following line of arguments to get a feeling for 
the expected T dependence of rj. (i) For V > T > Op, T F 
the electron fluid is a highly correlated classical fluid. 
There are many examples of such fluids - i ndeed , most 
classical fluids fall in this regime (jFrenkell . 11946 ) . The 
viscos it y of classical fluids i s wide ly observed (jFrenkel 
( 1946 ): iFerrer and Kivelsonl (|1999l )) to be an exponen- 
tially increasing function of decreasing temperature, (ii) 
For T F « T « J] P = \/EpV the fluid is still quantum 
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but not degenerate, and is still strongly corr elated. It 
has b een conje ctured on theoretical groun ds (| Andreev . 
19791 ) (see also ISpivak and Kivelsonl (120061 )1 that in this 
regime 



r,(T) cx 1/T. 



(13) 



III. THEORETICAL INTERPRETATIONS OF 
EXPERIMENT 

We now discuss some of the attempts that have been 
made to interpret the corpus of experimental observa- 
tion summarized in Section 1 in terms of the various the- 
oretical results outlined in Section 2. As is clear from 
Section 2, there is presently no well controlled theory 
that treats non-pcturbatively both the strong correlation 
and the disorder effects, so all such attempts involve an 
extrapolation of results from small r s to large r s , from 
zero disorder to finite disorder, or both. Thus, while we 
present arguments both in support of and against vari- 
ous proposed interpretations, none of our conclusions are 
irrefutable. 



A. Explanations based on classical "Drude" formulas 

A systematic attempt to explain a wide range 
of the experiments presented in Section 1 has 
b een undertaken b y 



Gold and Dolgopolov 



( 20001) , and 



Stern and Das Sarma 



(1985) 



19861). iDolgopolov and Gold 



Das Sarma and Hwang (1999 



2000 



2003, 



20041 ) employing the classical formulas for the resistivity 
which are valid for weak scattering and r s <C 1 and 
extrapolating them to the case r s 3> 1. Manifestly, 
this approach involves extrapolating results from the 
weak interaction regime into a regime in which they are 
strong. Since the electron impurity scattering is treated 
perturbatively (Born approximation) and the electron 
screening is computed at the level of RPA, there is 
no formal justification for the approach when r s > 1. 
However, the appeal of this approach is that it leads to 
explicit expressions for a wide range of physical quanti- 
ties which can be directly compared with experiment. A 
number of striking quantitatively successful comparisons 
be tween this theory and experiment have been rep orted 



bv lDas Sarma and Hwang) (|1998 



2000 



2003 1 120041) 



There are, however, several aspects of this extrapola- 
tion that we find troubling. 

The extrapolation to large r s does not simply involve 
quantitative shifts, but qualitative changes. Whereas at 
small r s , as discussed in Section flLAl both dp/dT < 
and dp/dBn < 0, the classical formulas extrapolated to 
large r s exhibit the opposite sign trends, dp/dT > 
and dp/dBn > 0. The change of sign at r s ~ 1 does, 
admittedly, bring the results into qualitative agreement 
with experiment on large r s systems, but it cannot be 
said to be a "featureless" extrapolation. Moreover, the 
origin of the sign change can be traced to the fact, 
already discussed above, that at large r s , the screen- 
ing length A sc obtained in RPA approximation is para- 
metrically smaller than the spacing between electrons, 
X sc y/Trn = (l/4)r^ 1 -C 1. Screening lengths less than the 
distance between electrons are clearly unphysical, and we 
worry that the same is true of other extrapolated results. 

This same approach has been extended to explain the 
decrease of p(T) as a function of increasing T at T > 
Tmax ~ Tp, where the electron gas is non-degenerate. 
Here, the resistivity is computed using the kinetic the- 
ory of a weakly inter acting electron gas scat t ering from 
Coulomb impurities (|Das Sarma and Hwand . Il999f ). In 
the framework of this approach the decrease of the re- 
sistance is associated with the increase of the velocity 
and decrease of the scattering rate as T increases. As 
a result p cx 1/T. Again, semiquantitative success has 
been reported in the comparison of this theory with ex- 
periment. Since this theory neglects all correlations in 
the electron gas or between scattering events, it is well 
justified when G> 1 and r s <C 1 for all T > Tp, and 
even for r s ^> 1 at temperatures T ^> V. However, the 
experiments in question deal with the temperature range 
Tp < T < V ~ r s Tp in samples with r s 3> 1; in this 
regime, the electron liquid is still strongly correlated, al- 
though not quantum mechanically coherent. 

Additionally, Drude theory does not incorporate an 
insulating state, much less a metal-insulator transition. 
It seems clear to us that the anomalous behavior of the 
metallic state at large r s are related to the appearance of 
a metal-insulator transition in the same devices. 
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B. Interaction corrections to the conductivity in the 
"ballistic regime" Lj <C I. 

For r s <C 1 and at high enough temperatures that 
Lt <C I but low enough that T < Tp, the contribu- 



tion (Zala et 



2001al la) to the resistivity of electrons 
scattering from the Friedel oscillations induced by impu- 
rities is given by Eq. El At small r s , where |Fo| <C 1, 
this expression gives the opposite sign of dp/dT than is 
seen in experiments on samples with r s 3> 1 and G> 1. 
However, when Eq. [6] is extrapolated to large enough 
r s , where plausibly Fq < —1/4, dp/dT changes sign, and 
p(T) becomes a linearly increasing function of T. (Where 
there is valley degeneracy, as in the case of Si MOSFET's 
where there are two va lleys, the same sign reversal oc- 
curs wh en Fq < — 1 / 15J Zala et all l2001al lbf) . ) The same 
theory |zala et all . l2001a|b ) predicts a decrease in the 
magnitude of the triplet part of Eq. [6] by a factor of 3 
when the electron gas is spin polarized by the applica- 
tion of B\\ > B* . The reason is that the singlet and 
L z = triplet parts of the two particle propagator are 
unaffected by £>|| , while the remaining two components 
of the triplet are suppressed by B\\ . 

While this approach also involves an extrapolation, in- 
volving a sign change at r s ~ 1, there is certainly nothing 
unphysical about a Fermi liquid parameter with a sub- 
stantial magnitude: F < —1/4. However, there are other 
aspects of this explanation of the experiments presented 
in Section 1 that we find problematic. 

Firstly, this is a theory of corrections to the conductiv- 
ity; even when the result is extrapolated to r s ~ 1 and 
T ~ Tp, this approach only makes sense if the correc- 
tion is small compared to the Drude conductivity itself, 
Go- Since the correction is a non-analytic function of T, 
it can be the dominant contribution to dp/dT for suf- 
ficiently small T <C E-p. However, it is something of a 
stretch to interpret the large fractional changes in the 
conductivity seen in experiments on samples with r s 3> 1 
and G 3> 1 in these terms. Over a range of temperatures 
and fields, T ~ Ep and/or Bn ~ B*, the analytic vari- 
ations of the Drude conductivity (which were the focus 
of the theory discussed in Section IIII. A|) should gener- 
ally be at least as large as these corrections. In more 
physical terms, it is hard to imagine that the scatter- 
ing of the electrons from the induced Friedel oscillations 



of the electron density can make a larger contribution 
to the scattering cross-section than the scattering from 
the impurity itself. Therefore, we think that this pro- 
cess (even when extrapolated to large r s , and taking the 
most optimistic viewpoint) cannot explain changes in the 
resistivity by more than a factor of two. However, in ex- 
periment, p is observed to increase (admittedly more or 
less linearly with T) by a factor of 4 in large r s p-GaAs 
heterostructures, and by a factor of 10 in Si MOSFETs. 

Secondly, although this theory does predict a suppres- 
sion of the T dependence of the resistivity by a parallel 
magnetic field, dp{T)/dT\ B ^ B , / dp{T) / dT\ B ^ =Q < 1, 
the measured ratio at low T i s significantly sma ller than 



that predicted by the theory (jZala et ali . l2001al fbl) 



In short, we conclude that, although the re s ults o f 



the perturbation theorv ( Stern and Das Sarma ( 1985 


): 


Gold and DoIeodoIov 


(19861): 


DolgODolov and Gold 


(2000 


); 


Zala et al. 


2001aib)) are likelv relevant to 



experiments performed on samples with relatively small 
r s and at relatively small values of B\\ and T, they 
cannot explain results of experiments at r s 3> 1 where 
almost all effects are of order one or larger. 

C. Scaling theories of the metal-insulator transition in 2D 

The scaling theories discussed in the previous section 
imply the existence of a metal-insulator transition at a 
sample dependent critical density, n = n c . At critical- 
ity, the resistance approaches a finite limit, p c , in the 
T — > limit which is probably universal. Strong n — n c 
and T dependences of the resistivity near criticality are 
governed by universal critical exponents and scaling func- 
tions. The appeal of such an approach is the robustness 
implied by universality. However, our principle interest 
is with the behavior of samples over a broad range of r s . 
Thus, other than to stress the importance of the exis- 
tence of the transition itself, we have chosen not to focus 
particularly on experiments very close to criticality. On 
the basis of the data in Fig. [TJ we identify an empiri- 
cal value for the critical resistance, p c f=a l(h/e 2 ). Thus, 
we can identify samples with low temperature resistances 
p > p c as "insulating" and with p < p c as "metallic" . 

Any theory which is based on electron interference ef- 
fects, in practice predicts substantial fractional changes 
of p(T, B») over plausibly accessible ranges of T only in 
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the near vicinity of the critical point. In particular, where 
G ^> 1, only fractionally small logarithmic variations are 
predicted. In contrast, in all systems which we describe 
in this review, the transport anomalies take place in a 
wide range of the electro n densities: up to five times the 
critical density (see, e.g. 



Altshuler et al 



(|200rt n. There- 



fore, most of the experimental data do not lie in the crit- 
ical region. The experiments deep in the metallic re- 
gion p(T, exhibit very large changes in the resistance 
(sometimes by the order of magnitude) even in samples 
where the conductance is as large as G ~ 10 — 20e 2 /h 
(see FigsHJU). I n f act > the temperature dependence of 
the resistivity significantly increases in samples which are 
farther away for the critical region. 

We therefore conclude that, although this theory quan- 
titatively describes experimental dat a on Si MOSFETs 



in th e close vicinity of the transition (jAnissimova et al 



20071 ) , it cannot explain large effects far from the transi- 
tion. 

D. Interpretations based on electronic microemulsions and 
the Pomeranchuk effect 

In an electron fluid with r s 3> 1, the interaction en- 
ergy is large compared to the kinetic energy, and so the 
short range correlations must certainly be Wigner crys- 
talline like, whatever the long-range emergent properties. 
No theoretically well controlled treatment of this prob- 
lem, capable of quantitative comparison with experiment, 
currently exists that treats on an equal footing the short- 
range crystalline correlations and the long-distance fluid 
character of the state. The nascent theory of electronic 
microemulsions is, however, an attempted first step in 
this direction. Here, we sketch the ways that many of the 
most significant qualitative aspects of the experiments 
can be understood from this perspective. 

1. Interpretation of the T and Bu dependences of the 
resistance for T < Tp 

The T and Bn dependences of the resistance at T < Tp 
can be qualitatively understood as consequences of the 
Pomeranchuk effect. As the fraction, /wc, of Wigner 
crystal grows with increasing T and Bu , this naturally 
produces an increasing p, i.e. a "metallic" T depen- 



dence and a positive magneto-resistance. Moreover, to 
the extent that the exchange energy, J <C fcrjT, /wc 
is a strongly non-analytic function which, for B\\ = 
and T <C Tp is linear in T . and for T <C B\\ is an ini- 
tially linearly increasing function of Bu which saturates 
at Bn > B* . Moreover, the T dependence of /wc is 
quenched by B\\ in the range T < By. 

As seen in Section 1, the experimentally measured p 
in samples with G 1 and r s 3> 1 exhibits all the same 
qualitative dependences as /wc, including the dramatic 
1-2 order of magnitude decrease in dp/dT produced by a 
field B \\ > B* . The t r anspo rt theory of microemulsions 
(jSpivak and Kivelsonl . [2006) which relates p to /wc is 
complex and incomplete. In the limit that the disorder 
potential is smooth and /wc 1, the dominant contri- 
bution to the resistance comes from electrons scattering 
from rare islands of WC, which are themselves pinned at 
minima of the disorder potential. In this case, it is easy 
to see that p oc /wc- 

For other forms of disorder, and in all cases where /wc 
is not small, the dependence of p on /wc is more com- 
plicated, although, presumably p is still a monotonically 
increasing function of /wc- Thus, the Pomeranchuk ef- 
fect gives a plausible explanation of the giant positive 
magneto-resistance in the "insulating" regime, p > h/e 2 , 
of the sort shown in Fig|T2J It also gives rise to the possi- 
bility that samples with p > h/e 2 can still have a metal- 
lic temperature dependence dp/dT > 0. Such behavior 
is occasionally seen, as shown in FigfTTl although so far 
only relatively small effects in samples with p no larger 
than ~ 3/i/e 2 . However, in principle, the Pomeranchuk 
effect can produce arbitrary large effects. In particular, 
on a qualitative level it can explain the existence of a 
metal-insulator transition as a function of Bu . 



2. Interpretation of the magnteto-resistance in a perpendicular 
magnetic field 

Moving to the effect of a perpendicular magnetic field, 
the experimentally observed behavior can be quite com- 
plex, as shown in Figsj6] [13] and [14] As is confirmed in 
tilted field experiments (also shown in the Fie fM]) . this 
complexity reflects a combination of spin and orbital ef- 
fects. The strong negative magneto-resistance seen at 
large fields, where the spins arc fully polarized, gener- 



24 



ally reflects the existence of quantum Hall phases, as dis- 
cussed in Section III. Gl For intermediate values of B± , 
where spin physics is important, the Pomcranchuk ef- 
fect can readily account for the existence of a large posi- 
tive magneto-resistance, as discussed above. (The small 
negative magneto-resistance, sometimes seen in metal- 
lic samples at small values of B± (see Fig|6]), is pre- 
sumably due to weak localization effects.) This gen- 
eral analysis applies either to samples in which the zero 
field resistivity is less than h/e 2 (Figf5|) or greater than 
h/e 2 (FigsfT51 and [T4]) . In the more resistive samples, 
a part of the negative magneto-resistance is presumably 
due to interference corrections to variable range hopping 



(jShklovskii and Spivakl . 



1991 



Zhao et al 



19911). 



3. Interpretation of the T-dependence of the resistance at 

T>T F 

The existence of a broad range of temperatures with 
Tpr s ~ V > T > Tp in which the electron fluid is still 
strongly correlated is one of the most clearly inescapable, 
and at the same time widely neglected features of the dis- 
cussion of the 2DEG with large r s . The best understand- 
ing we have of the transport theory, here, is obtained by 
arguing in analogy with other strongly correlated fluids. 
As discussed in Section III. H. 31 because there is no co- 
herence, the Fermi statistics of the electrons is relatively 
unimportant, but so long as Tp-^/r^ > T > Tp, quantum 
effects could be crucial. Moreover, because the interac- 
tions are strong and the samples are "clean," it is plau- 
sible that p oc 77, where 77 is the viscosity of the electron 
fluid in the absence of disorder. 

The best analogy is with the viscosity of liquid 3 He 
and 4 He in the same temperature range. The experi- 
mental fact t hat p ~ 1/T in t his temperature interval 



(see Figl2land lMills et ali (|l999T j) is not inconsistent with 
the conjecture, Eq. Q2J that the viscosity of liquid He in 
the same "se miquantum" re gime, obeys 77 ~ 1/T. (Mea- 



sured values (|Andreevl . ll979l ) of the viscosity in bulk 4 He 
are consistent with 77 ~ 1/T, but as far as we know, no 
comparable data exists in bulk 3 Hc, nor in 3 He or 4 He 
films.) 



E. Interpretation of experiments based on percolation 

Percolation is a classical concept whose relevance to 
quantum systems requires a separation of spatial scales 



(Efros and Shklovskii 



19841 ). Specifically, in order that 



different (randomly distributed) regions of the sample 
can be treated as sufficiently macroscopic to be character- 
ized by a local value of the conductivity, the correlation 
length of the scattering potential must be large compared 
to the relevant microscopic lengths needed to define a lo- 
cal electronic phase. At the very least, percolation is a 
useful concept only if the disorder potential varies on dis- 
tances large compared to the distance between electrons. 
In particular, since the quantum dephasing length, di- 
verges as T — > 0, percolation is never justified at very low 
T because it neglects the localization, interference effects, 
etc. in the "metallic" state. So percolation at best pro- 
vides a valid description at intermediate temperatures. 
However, if the disorder has sufficiently long-range cor- 
relations, percolation effects can readily mimic the finite 
T appearance of a metal-insulator transition. 

Such a theory can be developed in a form which is 
applicable to samples with small r s in the presence of a 
smooth scattering potential. (This situation can be real- 
ized, for example, in 2D samples where charged impuri- 
ties are separated from the 2D electron liquid by a wide 
spacer.) In this case, to describe the tr ansition one has to 



take into acc ount non-linear screening (jEfros et 
Foglerl . 120041 ). 



1993: 



Das Sarma et 



( 20051) 



Manfra et al 



( 20071 ) applied 



the same sort of theory, extrapolated to large r s , to ex- 
plain experimental data for low density electron liquids in 
in the metallic state near the point of the metal-insulator 
transition. Again, despite the reported qualitative suc- 
cesses of this theory, we believe that there are significant 
problems with its straightforward application in the cur- 
rent context. For example, as B» increases the Fermi 
energy of weakly interacting electron liquid increases, ul- 
timately by a factor of 2. Generally, this would be ex- 
pected to increase the faction of the sample that is in 
the metallic phase. Therefore such a theory would result 
in a resistance which is a decreasing function of increas- 
ing Bu , with even the possibility that it could induce an 
insulator to metal transition. Conversely, the experimen- 
tal data presented in Fig. 4 show a field-driven metal to 
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insulator transition. 



F. Experiments on double-layers 



Indeed, the assumption of very long-range correlations 
of the disorder is problematic in itself. Certainly, in Si- 
MOSFETs, it is generally accepted that the disorder has 
rather short-range correlations. More generally, the dis- 
tance to the remote dopants is rarely significantly longer 
than the typical distance between electrons. Thus, the 
large separation of length scales assumed in a simple per- 
colation analysis cannot be ubiquitous among the devices 
studied. 

On the other hand, the metal- Wigner crystal transition 
(and associated intermediate micro-emulsion phases) in 
the presence of weak disorder yields a state which can 
naturally be described by percolation through the metal- 
lic fraction of the system. Even for short-range cor- 
related disorder, the familiar disorder broadening of a 
first order transition can lead to "I mry-Ma" domains of 



large siz e (llmrv and Wortid (|1979h 



Aizenman and Wehr 



(|1989l ); Bakerl (|l993r )). Although, at the end of the 
day, a percolation picture based, in this way, on dis- 
order induced phase coexistence has many features 



I Das Sarma et al 



Manfra et al 



. 2005: 


Efros et al.. 


1993; 


Foglerl 


2004 



20071 ). above, it has the important differ- 
ences that: (i) The large length scales that justify the 
percolation analysis are self-organized, rather than in- 
herited directly from the disorder potential, and (ii) The 
resulting pattern of insulating and conducting phases is 
a strong function of any parameter, especially T and £?m , 
which affects the thermodynamic balance between the 
two phases. 

As we have already discussed, this general type of the- 
ory is in qualitative agreement with much of the exper- 
imental data. For example, due to the Pomeranchuk ef- 
fect, the fraction of the Wigner crystal increases with 
increasing B\\, which in turn can lead to a field-induced 
transition from a metal to an insulator. Another ad- 
vantage of consideration of this regime is that it offers 
a simple explanation of why highly conducting samples 
can still have such strongly T and H\\ dependent conduc- 
tivities. Unfortunately, however, no quantitative theory of 
the transport in this regime currently exists. 



We now tu rn to the anomalies in t he measured drag 
resistance po ijPillarisettv et aZ.l . l2003h in p-GaAs double 
layers with small resistances and large r s . 

A straightforward application of the Fermi liquid based 
theory valid at r s <C 1 yields results i n qualitative con 



tradiction with experimental results of iPillarisettv et al 



(|2003l ) . For example since in the framework of the Fermi 
liquid theory the electron scattering rate is inversely pro- 
portional to the Fermi energy, it natural to expect that 
the drag resistance decreases as B\\ increases. This is in 
contradiction with the experiment. 



It was shown bv lLevchenko and Kamenevl (|2008l ) that 
interference corrections to pd are increasingly impor- 
tant at low T and indeed do not vanish as T — > 0. 
This is in loose qualitative agreement with the fact that 
the experimental values of the drag resistance are much 
larger than the conventional theo retical values (jPricd . 



19831 : IZheng and MacDonaldl . 1993 ) , which vanish in pro- 
portion to T 2 as T — > 0. Secondly, according to 
Levchenko and Kamenevl (|2008l ). pd ~ p 3 - This is also 



in qualitative agreement with experiment where both 
Pd(B\\) and p 3 (B\\) increase with Bu and saturate at 
Bii > B*. On the other hand, as we have already dis- 
cussed, since this theory was derived in the limit r s -C 1 
and k-pl 3> 1, it cannot answer the more basic question of 
why p(T,B\\) itself increases as a function of T and Bu, 
nor why there is a maximum of p(T). 

A focused study based on an extrapolation of Fermi liq- 
uid based form ulas to the region where r s 3> 1 has been 



carried out by iDas Sarma and Hwan g (2005 ) in order 



to add ress the drag resistance data of 



Pillarisettv 



et al 



(|2003l ). According to lDas Sarma and Hwang! (|2005l ). the 
large enhancement of the magnitude of the drag resis- 
tance arises from a combination of a number of separate 
(somewhat technical) features of the specific experimen- 
tal system in question. Moreover, these authors find that 
both the resistance of individual layers, and the the drag 
magneto-resistance turn out to be increasing functions 
of Bu . As we have mentioned our concern here is with 
the underlying assumptions in this approach: Despite 
its many succes ses, the formulas for th e elect ron cross- 
section used by IDas Sarma and Hwang] <|2005f ) were de- 
rived under an assumption that the standard expressions 



2G 



for electron screening continue to apply when the screen- 
ing length, X sc is much smaller than the inter-electron 
distance. Specifically, in this approach, the change of the 
sign of the drag magneto-resistance which occurs on ex- 
trapolating the small r s results to large r s can be traced 
to the fact that X sc -C rT 1 / 2 , which is clearly unphysical. 

An attempt t o exp lain the experimental data 
(|Pillarisettv et all I2003T ) based on the assumed ex- 
istence of microemu l sion phases was carried out by 
Spivak and Kivelsonl (|2005f ). The main assumption is 
that a phase consisting of mobile bubbles of Wigner crys- 
tal embedded in the Fermi liquid is responsible for the 
large drag resistance. Since the bubbles have different 
electron densities than the surrounding Fermi liquid, they 
produce a large electric potential which is seen by the 
electron liquid in the second layer, thus producing a large 
Pd- However, the existence of mobile bubbles at the rel- 
evant temperatures is an assumption that has not yet 
been tested, since no reliable theoretical estimates have 
yet been made of the characteristic size of the bubbles, 
nor the range of densities over which the bubbles are 
stable, and no direct experimental imaging of an elec- 
tron microemulsion has yet been achieved. Once this 
assum ption is accepted, all qua litative features of exper- 
iment (jPillarisettv et all 120031 ) can be explained as con- 
sequences of the Pomcranchuk effect. 

G. n-dependence of B* 

In discussing the n dependence of B* , it is again 
useful to compare it to the analogous problem in 3 He. 
The problem of the density dependence of the satura- 
tion magnetic field and enhancement of the spin sus- 
ceptibility in the strongly correlated Fermi liquid 3 He 
near the crystallization tran sition has been discussed by 
Castang and Nozierea (| 19791) . There, it was pointed out 
that two different scenarios can be imagined with differ- 
ent consequences for the evolution of the magnetic re- 
sponse: 

(i) The system is nearly ferromagnetic, which means that 
the linear spin susceptibility is enhanced compared to its 
non-interacting value, but the saturation field, B* , is not 
suppressed. This would mean that the system is close to 
a Stoner instability. 

(ii) The system is nearly solid. In this case, both the lin- 



ear and all non-linear susceptibilities arc enhanced, and 
B* is suppressed compared to the non-interacting values. 

It seems to us that the data in the large r s 2DEG are 
more generally consistent with the second scenario. 

IV. CONCLUSION 

We have summarized a large body of experimental data 
and theoretical arguments which strongly imply the ex- 
istence of qualitatively new ( "non- Fermi liquid" ) physics 
of the 2DEG in the limit of large r s and weak disorder. 
These phenomena include, but are not limited to those 
associated with a metal insulator transition. We have 
also critically discussed some of the attempts to make 
contact between theory and experiment concerning these 
behaviors. 

The primary purpose of this article is to bring into fo- 
cus the problems of physics which are unresolved, and 
to stimulate further study, especially experimental study 
of these problems. Because of the absence of a well 
controlled theory of the metal insulator transition in a 
strongly correlated 2DEG, there has been a large amount 
of controversy surrounding the interpretation of these ex- 
periments. However, it is important to stress that there 
is no controversy concerning the experimental facts, espe- 
cially given that similar phenomena are seen by multiple 
groups in a variety of very different types of devices. 

We conclude our discussion with a brief list of some of 
the properties of the 2DEG which could readily be mea- 
sured, but which have not (to the best of our knowledge) 
been seriously studied to date. More generally, the widely 
celebrated progress that is ongoing in the fabrication and 
broader distribution of increasingly high mobility Si and 
(especially) GaAs based devices will open the possibil- 
ity of probing the intrinsic properties of clean 2DEGs 
at still lower densities and temperatures, where presum- 
ably all the previously observed behaviors will have still 
larger amplitudes. Moreover, as high mobility devices 
made with other semiconductors become available, fur- 
ther tests of the universality of the phenomena, as well 
as interesting particulars associated with differences in 
band-structure and the like, could add to our knowledge 
of this subject. 

(i) There is surprisingly little data concerning the prop- 
erties of strongly correlated electron liquids in the semi- 
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quantum regime Ep < T < ^/r^E-p ~ ^/EpV, where it 
is non-degenerate, but still highly quantum mechanical, 
and in the classical regime, ^fr^Ep < T < r s Ep ~ V, 
where it is still highly correlated but quantum effects 
are negligible. For example, experimental data on the 
magneto-resistance in a parallel magnetic field at T > Ep 
have not been reported for any of the materials in ques- 
tion. (It is also worth noting that systematic studies of 
the viscosity of 3 He, especially 3 He films, in the semi- 
quantum regime do not, to the best of our knowledge, 
exist and would be useful for comparison.) 

(ii) As the resistance is the most directly measured 
property of the 2DEG, most empirical information about 
its character comes from transport measurements. How- 
ever, thermodynamic propertie s such as compressibi l- 
ity can be measured on 2DEGs (jEisenstein et all Il992f ). 
There have been some recent measurements of com- 
pressibility showing a drastic difference between th e 



Allison et al 



metallic and i nsulat in g phases foultz and Jiangl (|2000t k 
IpOod) : lllani et all feoool )). It would be 



desirable, however, to perform more studies in sam- 
ples i_with_ev§n_j3igher used 



bv iDultz and Jiangl |200(j| ) 



Hani et al 



Allison et al 



|2006( 1. and 



( 20001 ) to elucidate the difference between 



high r s samples with strong metallic transpo rt and the 



2DEG's with lower r s (Eisenstein et 



19921) 



(iii) Thermoelectric and Nernst effects can reveal the 
strongly correlated nature of an electron system. How- 



of papers ( 


Molr 


ovan et al. 


Possanzini et al. 


( 


2004 


)) r< 



Faniel et al 



(2007): 



surements in strongly correlated 2DEGs, while the 
Nernst effect has never been measured in devices with 
large r s . 

(iv) Since drag effects are uniquely sensitive to charge 
density fluctuations in the 2DEG, further experiments in 
double layer systems, especially in the presence of per- 
pendicular, and (even more importantly) parallel mag- 
netic fields, could be very illuminating. For example the 
drag resistance in the semi-quantum regime T > Tp has 
never been measured. 

(v) While isolated experiments exist on the quantum 
Hall effect in devices with moderately large r s , systematic 
studies of the evolution from the quantum Hall states to 
the Wigncr crystal state as a function of increasing r s do 



not exist. Moreover, as discussed above, there is much 
to be learned about the phases and phase transitions at 
B±_ = by following the evolution of the various quantum 
Hall phases to small B± in clean devices with large r s . 
(vi) It would b e very desirable to extend the few ex ist- 



Zhu et al 



20031 ) on 



ing experiments (jLillv et ali . 12003c 
low density, large r s electron gases in ultra-clean n-GaAs 
devices, and in particular to measure effects of the mag- 
netic field. 

(vii) One very promising avenue for obtaining more lo- 
cal information about the nature of the 2DEG at large 
r s , more or less free of the complications due to quenched 
disorder, is to study their properties in mcsoscopic de- 
vices such as "point contacts" . Up until now, however, 
they have primarily been studied only in three regimes: 
the regime where electr on interactions are not important 



van Wees et al 



19881 ). the Coulomb blockade regime 



and the Kondo regime ( 


Gold 


haber-Gordon et al. 




1998; 


Kouwenhoven and Glazmanl. 


2001 


) in which there is a 



single localized electronic state through which tunnel- 
ing between the two metallic reservoirs occurs. How- 
ever, there is another interesting regime when there is 
a relatively large "depletion region" between two metal- 
lic reservoirs with high electron densities. In the deple- 
tion region, there is a strongly correlated electron liquid 
with low electron density whose density can be varied 
by changing a gate voltage. In this case as a function 
of gate voltage, aspects of the electronic micro-emulsion 
phases can be directly probed on a mesoscopic scale. In 
this context we would like to mention that, as discussed 



by ISpivak and Kivelsonl (|2006l ) , the Pomeranchuk effect 
can produce significant T and B\\ dependences to the re- 
sistivity through the depletion region. It is interesting 
that it can produce dependences of p(T,Bu) which can 
mimic some of the behavior traditionally associated with 
the Kondo effect. 

viii) There are va r ious new experimental methods 



Tessmer et al 



(1998): 



Hani et al 



(2001): 



Sciambi et al 



20081 )) being developed which hold the promise of pro- 
viding spatially resolved images of the evolving physics 
of the 2DEG at large r s . Needless to say, such data could 
revolutionize our understanding of the strong correlation 
effects in the 2DEG. 

(ix) Finally, although it may be a while until suffi- 
ciently well ordered materials are available for these pur- 
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poses, it is likely that advances in the study of double- 
layer graphene and related materials may open unprece- 
dented opportunities to study the properties of the 2DEG 
at large r s . Due its Dirac spectrum. r s does not depend 
on n in a single layer of graphene, and more generally cor- 
relation effects are quite different than in problems with a 
quadratic dispersion. However, bilayers of graphene have 
a quadratic dispersion, so r s — ► oo as n — > 0. 
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